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We generalize the quantum optical input-output theory developed for optical cavities to ultracold 
fermionic atoms confined in a trapping potential, which forms an "atom cavity". In order to account 
for the Pauli exclusion principle, quantum Langevin equations for all cavity modes are derived. The 
dissipative part of these multi-mode Langevin equations includes a coupling between cavity modes. 
We also derive a set of boundary conditions for the Fermi field that relate the output fields to 
the input fields and the field radiated by the cavity. Starting from a constant uniform current 
of fermions incident on one side of the cavity, we use the boundary conditions to calculate the 
occupation numbers and current density for the fermions that are reflected and transmitted by the 
cavity. 
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I. INTRODUCTION 



cavity in terms of the field incident on the cavity. 



In light of the remarkable achievement of Bose-Einstein 
condensation in 1995 |Q, there has been a growing ap- 
plication of ideas from quantum optics to matter waves. 
This new field of atom optics Q has included both the- 
oretical and experimental investigations of matter wave 
coherence (|, |, [| ||, Q], atom lasers J|, ||, non-linear 
effects in matter waves including matter wave mixing 



jl0|, 11 , parametric amplification and squeezing in cou- 
pled optical and matter waves Jl^j . However the exten- 
sion of these ideas to degenerate Fermi gases has proven 
difficult because of the Pauli exclusion principle, which 
prohibits one from developing simple theoretical models 
based on only a few normal modes of the Schrodinger 
field. In addition to this, Fermi fields do not possess a 
classical limit analogous to the coherent state for Bose 
fields thereby making it impossible to develop semiclassi- 
cal mean-field theories such as the Gross-Pitaevskii equa- 
tion for Bose fields. 

All in all, the physical intuition obtained from quan- 
tum optics cannot be directly applied to theoretical in- 
vestigations of fermions. It therefore seems necessary 
that in order to make progress in the theory of fermionic 
atom optics, fundamental model systems in quantum op- 
tics need to be reanalyzed from first principles. Recent 
work in this direction indicates that four-wave mixing 
and coherent amplification of matter waves can occur 
in fermionic systems as a result of cooperative many- 
particle quantum interference analogous to Dicke super- 
radiance [H [jj. 

The purpose of this paper is to consider another model 
system, the atomic analog of an optical cavity with two 
partially transmittive mirrors. A schematic of our sys- 
tem is illustrated in Fig. 1. It consists of an atom cavity 
formed by two potential barriers with a finite number of 
bound states. The cavity states are coupled to a contin- 
uum of free particle states on either side of the cavity via 
tunnelling through the barriers. Our goal is to develop 
an input-output theory for fermions in the atom cavity 
that allows one to calculate the field radiated out of the 



The input-output theory for a single mode of a lossy 
optical cavity was developed by Collett and Gardiner in 
the form of quantum Langevin equations for the cavity 
mode |[q |. The great utility of this theory is that it al- 
lows one to incorporate the effects of quantum noise on 
the output field transmitted by the cavity as well as in 
the intracavity dynamics. Collett and Gardiner's formal- 
ism has been extended to bosonic matter fields in order to 
model the output coupling of atoms from a Bose-Einstein 
condensate in a single mode of an atom trap . As we 
will show below, the necessity of treating all modes of 
the atom cavity for fermions leads to novel features not 
present in the single mode bosonic theories. Most sig- 
nificantly, the eigenstates of the cavity become coupled 
due to their mutual interaction with the same external 
continuum states. Secondly, the coupling of the reser- 
voir modes to all cavity modes leads to the creation of 
coherences between fermions occupying different single 
particle modes in the radiated field even if the incident 
field is completely incoherent. 

The outline of the paper is as follows: In section II, we 
present our physical model for an atom cavity coupled to 
a continuum of reservoir states. In section III, we derive 
a set of quantum Langevin equations for the fermionic 
annihilation and creation operators of the eigenstates of 
the cavity in terms of both the input and output fields. 
These results are generalized to a two sided cavity in 
section IV. In section V, we consider a constant current 
of fermions incident on one side of the atom cavity and 
calculate the steady state statistics of the fermions trans- 
mitted through the other side of the cavity. In appendix 
A, we show that in a manner similar to the bosonic case, 
the presence of noise operators in the Langevin equa- 
tions is necessary to preserve the anticommutation rela- 
tions for the fermion operators. In appendix B we derive 
explicit forms for the coupling constants, which connect 
the intracavity modes to external continuum modes via 
a tunnelling Hamiltonian. 
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FIG. 1: Schematic diagram of ID atom cavity system. 



II. PHYSICAL MODEL 



A physical schematic of our system is illustrated in 
Fig. 1. For simplicity, we restrict ourselves to one spatial 
dimension, which allows space to be divided up into five 
distinct regions. The region —a<x<a between the two 
potential barriers of height Vq represents the atom cavity. 
For thick barriers, the number of bound states of t he 
cavity, N + 1, is given by Avr/2 < = yj2mV a 2 /h 2 < 
(N+1)tt/2 where m is the mass of the atoms ]l8| . We will 
focus on the case where N 3> 1. The regions —L — a — d < 
x < —a — d and a + d<x<L + a + d represent the 
left and right reservoirs, respectively. L is the length 
of the reservoir region. We let it go to infinity so that 
the fermions are described in terms of a continuum of 
free particle plane wave states. Atoms located in the 
reservoir regions with energies less the Vq couple to the 
cavity states by tunnelling through the potential barriers 
located at —a — d < x < —a and a < x < a + d. 

Since the wave functions for atoms with energies 
greater than Vq are not spatially localized in either the 
cavity or the reservoir regions, we restrict ourselves to 
single particle states with energies less than Vo- In this 
case we can meaningfully speak of left/right reservoir 
states and cavity states since the single particle wave 
functions decay exponentially inside the potential bar- 
riers. Such a restriction is valid provided the initial 
state does not contain any occupied states with energies 
greater than Vq and two-body collisional interactions be- 
tween atoms, which can cause atoms to be scattered into 
higher energy states, are negligible. The latter condi- 
tion will indeed be satisfied for ultra-cold spin polarized 
fermions since s-wave collisions are forbidden and p-wave 
collisions are negligible at these temperatures. Under 
these conditions, the states with energies larger than Vq 
are not coupled to states with energies less than Vo- 

The second quantized Hamiltonian for the cavity- 
reservoir system in the subspace of states with energies 



where H$, Hl, and Hr are the free Hamiltonians for the 
system (i. e. the atom cavity) and the left and right 
reservoir states, respectively, 



N 



H S = ^ h^ncicn 
ri=0 

H L = hw k a\a k 

k 



k 



hu} k b\b k - 



(2) 
(3) 
(4) 



Here, c„ is a fcrmionic annihilation operator that de- 
stroys an atom in the cavity with wave function <f>n (x) 
and energy M7„ = h 2 K 2 /2m. Similarly, a k and b k are 
fcrmionic annihilation operators that destroy an atom in 
the left and right reservoirs, respectively, with the wave 
function 

4> k l ' r \x) =exp(ik{x±(a + d))/L 1/2 

and energy hio k = h 2 k 2 /2m in the regions outside the 
barrier. 

The coupling between the system and the reservoirs, 
Hsr and Hsl, is given by effective tunnelling Hamilto- 
nians [fil Wl m m 
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In all cases the summation is restricted to those states 
with energies below the barrier. Explicit expressions for 
the tunnelling matrix elements, n nyk and K n)k , are given 
in Appendix B. We note here that in 1-D, the coupling 
constants depend only on the magnitude of k and not its 
sign. 

In contrast to quantum optical systems, which are of- 
ten approximated as a single cavity mode with a large 
occupation number, a full multi-mode treatment is re- 
quired for fermions even if the number of fermions in the 
cavity is small (~ 1). This is because of the Pauli prin- 
ciple, which forbids more than one atom from occupying 
the same cavity state, and thereby prevents one from sin- 
gling out a particular state as being more important than 
the rest. 

We conclude this section by noting that the general re- 
sults presented below do not depend on the precise nature 
of our physical model. We use a stepwise constant po- 
tential because it leads to simple analytic results for the 
coupling between the reservoirs and cavity states. Our 
model system in Eq. (|l|) could be applied to any fermion 
system in which a finite number of discrete states are 
linearly coupled to a dense continuum of states. 
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III. SINGLE SIDED CAVITY 

We now derive a set of integro-differential equations of 
motion for the cavity operators that only involve the ini- 
tial or final state of the reservoir operators. We proceed 
by formally integrating the equations of motion for the 
reservoir operators and substituting these back into the 
equations for the cavity operators. In this section, we set 
Hji = Hsr = so that the atom cavity is coupled to a 
single reservoir. This is analogous to an optical cavity 
with a single partially transmittive mirror. The inclu- 
sion of the right reservoir is summarized in the following 
section. 

It is convenient to work with slowly varying operators 



in the interaction representation, 

c„(i) = e-^c^t), o fc (t) = e^a k (t). (7) 

By formally integrating the Heisenberg equations of mo- 
tion for a k (t) from the initial time to to t, 

a k (t) = a k (t ) - <n,k f dt'e^^'c^t') (8) 

where a k (to) are the operators for the input field incident 
on the cavity barrier, and substituting this solution into 
the equations of motion for c n (t), we obtain 



dt 
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Here, a Il!TO (r) is a reservoir correlation function, 
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that decays to zero in a characteristic time r c due to the 
destructive interference between the different oscillations. 
Note that r c depends, in general, on the cavity states n 
and m that are coupled by a nj m(T), and r^ 1 is on the 
order of the bandwidth of the reservoir, Vq/H. Further- 
more, if we assume that c m (t) only changes significantly 
over a time scale T m 3> t c , then we can make the Markov 
approximation by setting c m (t — r) = c rn (t) in Eq. (|^). 
For time intervals t — to 3> t c , we can then make the 
replacement 
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These expressions are defined in the continuum limit that 
Ek^(L/2n)Jdk. 
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FIG. 2: Plot of a„,„(r) for n = 49 and d/a = 0.0001(solid 
line), 0.001(dotted line), 0.01 (dashed line). Times are mea- 
sured in units of Qg . Vo and a were chosen so that the cavity 
contains 50 bound states. 



The Markov approximation assumes that the correla- 
tion function decays very rapidly, which requires that 
K n,k vary slowly with k. Fig. 2, shows the correlation 
function |a n! „(T)| for the highest energy bound state in 
the potential well. This state has the longest correlation 
time since the | «n.,fc | 2 are the largest for this state and 
because the summation over reservoir states is restricted 
to states with energies less than Vo. For the three cases 
plotted in Fig. 2, the smallest value of among all 



the cavity states is O^ft^ 1 { d /a = 0.0001), O^ft^ 1 
(d/a = 0.001), and 1.36S1q 1 [d/a = 0.01). In each case 
we see that the correlation function goes to zero in a 
time much shorter than T n J„n- This shows that the 
Markov approximation is a very good approximation for 
our system with T n /rc > 10 2 . 

By combining the above results, we obtain a quantum 
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Langevin equation for each of the cavity modes, 

Cn(t) = -in n c n (t) - ^(7n,m + iAn,m)Cm(t) + F% n \t). 
m 

(15) 



In Eq. ([L5D, we have defined the input noise operator 
Fi m \t) as 



■UO h (t— to) 



a k 



(to) = y^^n.fcQfc ft) 



The operators Sfcfti) represent the modes of the output 
field that contain the field radiated by the cavity at earlier 
times. By substituting this expression into the equations 
of motion for c n ft), making the Markov approximation in 
the integrand, and transforming back to the Heisenberg 
representation one obtains 

C n (t) = -m n Cn(t) + ^(7n,m " lA„, m )c m ft) + F ( n ° ut) {t). 

m 

(21) 

(16) Here, Fi° ut) (t) is the output field noise operator for state 



where 



a fc m) ft) = Gfcfto) exp(— iwfeft - t )) 



is the annihilation operator for mode (j)^ (x) of the input 
Fermi field at time t, and the total input field operator 
is therefore 



That is, ^i tn '(x, t) is the free Fermi field that propagates 
from the initial time to to t in the Heisenberg picture. 

Since the initial reservoir operators obey the an- 
ticommutation relations {afcfto), at/ fto)} — $k,k' an d 
{afefto), dfc'fto)} = 0, it is easy to show that the noise 
operators obey the anticommutation relations, 



{F^(t),Ft ) Ht-r)} = e 



(17) 



{FtHt),F^\t-r)} = (18) 

Furthermore, if we restrict ourselves to time scales much 
longer than the correlation time r c , then we can approx- 
imate the correlation function in Eq. ( p"7| ) by a delta 
function times the area under a n ^ m (r), so that 



{Ft\t),F^\t~r)}^2 ln . m 5{r). 



n(m)t 



(19) 



Before proceeding, there are several features of Eq. 
( |i"5| ) that are worth mentioning. First, the dissipative 
term J2 m (ln,m + «A„, m )c m ft) gives a damping term plus 
an energy shift for n = m while for n ^= m there is a 
non-zero coupling between cavity states. The coupling 
between cavity states is a result of all states coupling to 
the same reservoir, which leads to an indirect coupling 
between cavity states. Second, the noise operators cou- 
ple all of the reservoir states to each of the cavity states. 
Moreover, the noise operators are different for each cavity 
state due to the n-dependence of the coupling constants. 

Instead of solving for the reservoir operators in terms 
of the initial time to, one can instead solve for a^ft) in 
terms of a final time t\ > t, 



fife ft) = fifcfti 



(20) 



F£»')ft) = 5> n>fc e- 



-ia> fc (t-ti) 



Ofeftl) = «n,fcOfc H*) 



(22) 



where we have defined the output field annihilation op- 
erator a^ wt \t) for the mode tf>^(x). The ft) are 
related to the output field of the reservoir by, 

*<»*>(M) = £«L Mt W0»O- 



It is clear that ^( out \x,t) represents the free Fermi field 
for the reservoir that propagates from t to the final time 
t\. It is easy to see that if the cifcfti) obey normal 
fermionic anticommutation relations, then the anticom- 
mutators for the output noise operators are the same as 
Eqs. dHH). 

The boundary condition for the barrier separating the 
cavity from the reservoir is obtained by subtracting Eq. 
(H) from Eq. @, 



Ft 1 if) 



fM(t)=2^ 7n , mCm (f) (23) 



Equation ([23|) relates the noise operator for the output 
field to the input noise operator reflected by the barrier 
and the field radiated by the cavity. It is of the same 
form as the boundary condition for an optical cavity |l6|| 
except that in our case there are separate noise operators 
for each cavity modes since we cannot, in general, assume 
that the coupling constants are independent of the cav- 
ity state. In Appendix A, we use Eq. ( ^3| ) and (15) to 
derive the anticommutation relations between the cavity 
operators and the noise operators at arbitrary times. 
Equation (E3T) is not very useful since it is the mode 



operators of the output field, *■* ft) , that are needed 
to calculate properties of the output field such as mode 
occupation statistics, current density, etc. Therefore, we 
must extract from Eq. (^3|) a boundary condition for 
the annihilation operators of the modes of the input and 
output fields. 

First we note that in the limit of an infinite potential 
barrier the system-reservoir coupling vanishes, K n> k = 0. 
In this limit, a fermion incident from the left is perfectly 
reflected by the barrier. Since ^ m \x,t) and ty(° ut \x,t) 



5 



are the free Schrodinger fields that propagate forward 
in time from t — > — oo to t and from t to ti — ► oo, 
respectively, it follows that 

a k (t )e^ to =-a_ fe (t 1 )e iu * tl . (24) 

The total field is the incident plus reflected fields, 

*(.t, t) = ^ {in) (x, t) + ¥ out ^ 0, t). (25) 

ft follows from @ and @ that *(-a - <f,i) = and 
that the eigenmodes of the reservoir are standing waves 
for Vq — ► oo. Note that ( p4| ) and (^5|) are second quan- 
tized versions of the relations for the incident and re- 
flected wave functions from an infinite potential barrier 
& 

For the reservoir, the displacement operator is given 

by 

D(x) = e- lPx ' h 

where P =-Y^k fi>ka k a,k is the momentum operator for the 
reservoir J23| | . Multiplying (^3|) by D (x) on the left and 
D^(x) on the right gives after some manipulation 

k 

= 2^ Kll ,,<^(l] m - Wt )cmWe fc . (26) 

m,k 

Multiplying Eq. ( ^6| ) by ^_^ a _ d dxe~ %k ' x and taking 
the limit L — > oo gives the boundary condition for the 
modes of the input and output fields, 

4 n \t) - a { °f{t) = 2nJ2 K * m<k d(n m - u k )c m {t) (27) 

m 

Physically, Eq. (|27]) says that the difference between the 
mode of the output field propagating away from the bar- 
rier with momentum — hk is the input field with momen- 
tum Hk reflected by the barrier plus the field radiated 
by the cavity. Furthermore, only those states of the cav- 
ity that have the same energy as the reservoir mode can 
radiate into that mode. For a one dimensional system, 

I 



One may also derive Langevin equations analogous to 
Eq. ( ^jf|) involving the output noise operators for the two 
reservoirs, which can be used along with (|34|) to derive 
the boundary conditions for the noise operators in the 



there will only be a single cavity mode that contributes 
to the right hand side of Eq. (|27|). The <5(f2 m — uik) comes 
from Eq. ( [Ti| ) where we assumed times much longer than 
width of the reservoir correlation function. Hence, it 
follows from the uncertainty relation AEAt ~ H that 
the range of reservoir energies that couple to each cavity 
mode goes to zero as t — t$ — * °o- 



IV. TWO SIDED CAVITY 

The generalization of the preceding to a two sided cav- 
ity, HsR-, Hn ^ 0, is straightforward since the reservoirs 
couple independently to the cavity. For the right reser- 
voir, we define the input and output noise operators, 

k 

Gl°" t) («)=E^ 6 i OUt) W ( 29 ) 
fc 

where 

6^ n) (t) = 6 fc (*b)e- <w * (t - to) (30) 

b ( r t \t) = b k (t 1 )e- iu ^ t - t ^ (31) 

represent the input and output annihilation operators for 
the free field modes of the right reservoir at time t. Asso- 
ciated with the right reservoir noise operators are damp- 
ing constants 

In.m = «n,fc«m,fc^(^m _ ( 32 ) 

k 

and radiative energy shifts 

A n m = P y K n ,fcK m k — . (33) 

t-f ' U m -u>k 

k 

The quantum Langevin equations for the cavity mode 
operators expressed in terms of the input fields from the 
left and right are then 



(34) 

I 

left and right reservoirs, 

Gt\t) - G(° ut \t) =2j2%mc m (t) (35) 

m 

F (in) {t) _ F (out )(t) = 2 £ 7nimCm (t). (36) 



c n (t) = -in n c n (t) - E ((7n,m + i^n,m)c m {t) + (%, m + iA„, m )c m (t)) + F% n) (t) + G^ n) (t) 

m 



6 



Using the boundary condition b}Ato)e ZUJkto = 
—b-k{ti)e luJktl that corresponds to ( K4j) for the free 
fields in the right reservoir, the boundary condition for 



where 



J™) 

J k 

as Eq. 



.(out) 
J k 

One finds 



br>(t) and b'° ul \t) can be derived in the same manner 



b ( r\t) - b { °t t] {t) = 2^^n, k m m - w fc )cm(t). (37) 

m 

4 n \t) - a { °f(t) = 2tt£ K* m , k 5(Q m u; k )c m (t). (38) 



Equations (34), (J37|), and ( p8[ ) are the central result of 
this section. Given an initial state for the two reservoirs 
at to, Eq. ([mJ) can be used to calculate the state of the 
cavity at some later time. The mode operators for the 
output field of the left and right reservoirs can then be 
determined from the boundary conditions, (|37|) and (38). 



V. OUTPUT FIELD STATISTICS 

We illustrate how to utilize these results for a partic- 
ular initial state of the cavity plus reservoirs. Specifi- 
cally, we assume that the atom cavity initially contains 
no atoms and that the right reservoir is likewise in the 
vacuum state. Furthermore, the left reservoir contains a 
"beam" of fermions incident on the barrier at x = — a — d, 
with a spatially uniform current density equal to phq/m 
where p — M/L is the linear atomic density and Af is the 
total number of fermions. This physical configuration is 
represented by the initial state vector, 



!*(*>)) 



n 4(*o)io) 



(39) 



\k-q\<k f 



where kp — 2irp is the Fermi momentum. |^(to)) rep- 
resents a zero temperature Fermi distribution that has 
been given a Galilean boost that displaces the gas by q 
in k-space. This is analogous to the optical case where an 
incoherent white light source is used to drive an optical 
cavity. 

Even though | (to)) is a state with a fixed number 
of atoms, it acts like a constant input flux of fermions 
on the cavity. This is because of the implicit use of the 
Born approximation in the derivation of the Langevin 
equations, i. e. the reservoir is assumed to be so large 
that the back-action of the system is negligible. 

Using the results of the previous section we can cal- 
culate how the state of the left and right reservoirs are 
modified due to their coupling to the cavity. In particu- 
lar, we focus on the occupation numbers 



n! k R \t)=b k out) \t)b k out \t) 



(40) 
(41) 



as well as the current density operators for the reservoirs 



2mL 



2mL 



^(2fc + g)4 rat)t (t)4 + ^(t) (42) 

k 

^2(2k + q)b k ° u ^(t)b k °^(t) (43) 



are the spatial Fourier components of the current. Note 
that mL(jQ L,R \t)) is the average momentum in the out- 
put fields. 

Equation (|3~i| ) for c n (t) can be numerically integrated 
but we note that because of the exponential dependence 
of the reservoir-cavity coupling constants on the barrier 
height and thickness, the off-diagonal coupling is much 
smaller than the energy difference between the cavity 
modes, 



I a 



| --^ ~t~ ^fn,m \ : \ An,m ~t~ ■ 



for n m. In fact, a numerical evaluation of |7n,m+7n,m| 
and |A ni?71 + A„ jm | using the coupling constants of Ap- 
pendix B indicate that these terms are at least two or- 
ders of magnitude smaller than the energies of the cavity 
modes. It is therefore an excellent approximation to ne- 
glect all off-diagonal terms in the equations of motion. 
Furthermore, we can perform a renormalization of the 
cavity mode energies by absorbing the radiative energy 
shifts into them, O m + A mjlTl + A m ^ m — > £l m . 

For times much longer than the lifetimes of the cavity 
modes, (t — to)(7n,n + In.n) 3> 1, the system reaches a 
steady state with the solution, 



K m,k a k ^ (t) 

fe 



K 



7n,kUfc 



i(fi m -uj k ) + r.,. 



(44) 



where T m = ^ m ,m+lm,m- It is easy to see from Eq. ( |44| ) 
that the occupation numbers for the cavity modes have 
a Lorentzian profiles, 



\K m ,k\ 2 ( n k(to)) 



(45) 



where 



<n fc (* )> = e(k F -\k-q\) 



are the occupation numbers for the input field. Due to 
the broadband nature of the input field all cavity modes 
with energies in the range LUk F + q — u)k F - q will have a 
significant population with higher energy cavity states 
having larger populations due to the |/v m ,fc|'s exponential 
dependence on energy. 

Using Eqs. (|44|) and the boundary conditions |37]) and 
(|38|), we obtain steady state expectation values of the 
occupation of the output field modes , 
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r m A'(o m -^fc)|K ttl , fc |' 
(n m -^ fc ) 2 +r^ 



(n k {h))+^Y, r . 



2 5(fi m -!i„)J(!) m -u t ) 



(o m - Wfc ,) 2 +r^ 



»/f(*o)) (46) 



,(*)\ - 



l«m,fe| 2 |t m fc H 2 5(o m -n„)5(n m - Wfc ) 



(Or, 



2 +r2 



rife/ (t )). 



(47) 



Equations (46) and ( |47| ) represent the changes in oc- 
cupation numbers due to reflection and transmission 
through the cavity. The most significant feature of Eq. 
( [46| ) is that fermions in the input beam are perfectly re- 
flected by the barrier, (n^) = (nk(to)), unless there 
exists a cavity mode that is degenerate in energy with 
state k of the reservoir. The second term represents the 
interference term between the fermions reflected by the 
barrier and the fermions that have tunnelled into the 
barrier and subsequently tunnelled back out into the left 
reservoir. The last term in Eq. ( f46| ) represents the tun- 
nelling of atoms into mode — k as a result of atoms from 
mode k' that have tunnelled into the cavity and then 
tunnel out of the cavity. 

In order to gain additional physical insight, we simplify 
these expressions by noting that the denominator in Eq. 
(fl7]) as well as the last term in Eq. (|4^) are sharply 
peaked around u>k> = O m = u>k- Therefore we can replace 
k' with k in this term and drop the summation over kl . 



Using r m = 27„ 



27r|«; m ,fc| S(Q m — uok) we have 



l -ki 



where 



(n*(t )) - C(u k ) ((n fc (to)) - (n_ fc (to)» , (48) 



,(R)\ 



£(w fc )((n fc (to)) + (n_ fc (to))), (49) 



cupying the states k — 2tt/ L, (2tt/L) x 7501. Note 
that we have taken the reservoir to consist of discrete 
k states, k = 2 r Kn/L 1 with n = 0, ±1, ±n max and 
n m ax = L\/2mVo/2irh = 10 4 . Vq and a were chosen 
so that the cavity contained 50 bound states. Each line 
in Fig. 3 corresponds to a single reservoir k state and 
as such, the width of the lines are greatly exaggerated. 
The plots give good qualitative agreement with the above 
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FIG. 3: Plot of (a) (nj^(t)) and (b) {n ( k R) (t)) for d/a = 10 
and a/L = 1.2 x 10~ 5 . k is in units of 2-k/L. 



C(u) k ) 



L m 

(n m -cj k )2+T* 



rn 



(50) 



When (n-k(to)) = 0, Eq. (48) indicates that 



»<2) 



0, 

a result of the complete destructive interference between 
the fermions that are directly reflected by the barrier and 
the fermions that tunnel out of the cavity into the left 
reservoir. At the same time, Eq. ( f49| ) indicates that 
fermions resonant with a cavity mode tunnel through to 



the right side with unit probability, (n 



(R)\ 

±k j 



(nk(to))- 



These transmission resonances are similar to the situa- 
tion in an optical Fabry-Perot cavity. 

Fig. 3 shows a plot of (n^) and (n^) using Eqs. 



(Ha) and (147]) for an incident beam with fermions oc- 



discussion with each line located at the fc state that is 
closest in energy to a particular cavity state. The ampli- 
tude of the transmission resonances is close to 1 for the 
higher energy k states, while for the lowest energy res- 
onances, the amplitudes are about 0.25. The reduction 
in the amplitude for the low energy states in comparison 
to ( f48| ) and ( |49|) comes from T m = 2-K\n m ^\ 2 r)(p, m ) > 
27r|K mj fe| 2 (5(f2 m — u>k) where T}(u)) ~ lu^ 1 ^ 2 is the contin- 
uum density of states in the reservoir (see Eq. (|l 

The steady state current density in the output field of 
the left and right reservoirs are given by, 
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m 



(j-q ) = — hfi - }^( 2k + 1) I ~ 2lT 1^ K m,k+q K m,k ' 

k \ m 



and 



, , 2 V" l^m,kK* k+ gK* k ,K ntk 'S(Q m -U k )5(Q n -UJ k+ q) , . 

+4n 1^ — ( To — x , r w-/ n rX~Fl — («k'(*o)) -(51) 

^— ' {-i{ii m - ujk>) + r m )(t{u n - (Jv) + r n ) i 

Tilt - Tt ■ F\> 



4 > = ^—F 2J^ k + 1 4?r 2^ — / vo \ I -p uTTn rX~Fl — ( n k'(h)) , (52) 



respectively. 

The first term on the rhs of Eq. ( pl| ) is the inci- 
dent current reflected by the barrier. The average mo- 
menta in the output fields are proportional to (Jq ) — 

-H/mL J2 k fc(raLfe) and {jo ) = n / mL T,k k { n k )• In 
the left reservoir, (j^) ~ —phq/m since most of the 
fermions are perfectly reflected by the barrier. In the 
right reservoir, the transmitted current {jjj ) ~ since 
for those states that are resonant with the cavity (n k } = 

{n_u) while for non-resonant reservoir states, {n k ) ~ 0. 
Physically, this is due to the fact that atoms tunnel into 
the right reservoir from a standing wave cavity mode 
and therefore, they have equal probability to tunnel into 
states with positive and negative momentum. 

The terms are the spatial modulations in the 

current density that build up in the reservoirs as a result 
of the reservoir-cavity mode coupling. This can gener- 
ate a coherence between the k and k + q modes when 
there is a finite amplitude for an atom initially in state 
k to tunnel into the cavity and then tunnel back out of 
the cavity into state k + q. The change in (j^) is of 
order k 2 rather than k 4 as was the case for Eq. ( [46| ) 
since the current only involves generating a coherence 
between k and k + q rather than the transfer of popula- 
tion. Furthermore, the n 2 terms in ( |5l| ) are only finite 
for |wfc+g — uj k \ < r m , which implies that the coherence 
is only generated between reservoir states whose energies 
lie within the linewidth of a particular cavity mode. Con- 
sequently, decreasing the thickness and height of the bar- 
riers will make the lincwidths of the cavity states larger 
and thereby increase the magnitude of q ^ components 
of the current. 

Fig. 4 shows a plot of (j^) for several values of q. The 
q =/= components of the current are about two orders of 
magnitude smaller than {j^} and they decay away with 
increasing q. We do not plot (j^ ') since the current is 
equal to to within our numerical accuracy. This is be- 
cause the cavity linewidths are so narrow that it becomes 
nearly impossible to satisfy both the delta functions in 
the numerator and the Lorentzian denominators of Eq. 



0.03 




1 2 3 4_5 6 7 8 9 

q 

FIG. 4: Plot of (jjj } for the same parameters as Fig. 3. q 
are in units of 2n/L and the current is in units of h/2ma 2 . 
phq/m = 5.27 in these units. 

(H) for q ^ . 

VI. CONCLUSION 

In this paper we have extended the quantum op- 
tical input-output theory to atom cavities containing 
fermions. This formalism can easily be applied to intra- 
cavity nonlinear atom optical processes such as four-wave 
mixing between fermions |L3| or coherent photoassoci- 
ation of fermions into molecular bosons pa ]. In a fu- 
ture work we plan to extend this work to treat non- 
Markovian dynamics [^7| as well as including the effects 
of the spatially delocalized single particle state with en- 
ergies greater than Vq. 
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VII. APPENDIX A: ANTICOMMUTATION 
RELATIONS BETWEEN CAVITY AND NOISE 
OPERATORS 

In this appendix, we examine some of the consequences 
of Eqs. ([l5) and (p3|). Correlations between the input 
and output fields and the cavity modes may be expressed 
in terms of the anticommutators for the cavity and noise 
operators, i. e. nonvanishing equal time anticommuta- 
tors imply that the operators are not independent. From 
Eq. (|15|), one sees that the solution for c n (t) can be 
expressed in terms of the initial conditions for the cav- 
ity operators, c n (to), and the ak(to). It follows imme- 
diately that {c n (t),F^ n) (t')} = for all t and t' since 
the only non-vanishing anti-commutators at to are be- 
tween creation and annihilation operators. In a simi- 
lar manner, it immediately follows from Eq. (^l|) that 
{c n (t),F^ ut \t')} = 0. Formally integrating Eq. (§|) 
from to to t shows that c n (t) depends on (t') for 
t 1 < t (note that c n (t) will, in general, depend on all of 

the noise operators, Fm n \t'), not just n — m, due to 
the intra-cavity mode coupling). Using Eq. ( |l^ ) it then 
follows that 



{ C „(i-T),Fjr>t(t)} = o 



(53) 



for r > 0. This is nothing more than a statement of 
causality, the system can only depend on the past input. 
In a similar manner, formal integration of Eq. ( pl| ) from 
t + r to t\ implies that 



{c n (t + T),Fi°^(t)} = 



(54) 



for r > 0. Again, this is nothing more than the statement 
that the output of the system can only depend on the 
state of the system in the past. Using Eqs. (^34]) and 
(p3|), one obtains 

{c n {t + r), Ft )f (*)} = 2 ]T l* m , p {c n (t + r), 4(t)} 



while using Eq. (|53j) along with (J23J) gives, 
{c n (t-r),F^(t)} = -2V 7 ; p {cn(t-r),ct(f)} 



for r > 0. It follows from Eq. (|23|) that the equal time an- 
ticommutators can be written as {c n (t), Fm* 1 ^ (t)} — 
± J2 P Jm,p{ c n(t),4,(t)}+A n , m . The operator A n<m is de- 
termined from the equations of motion, ( |l5| ) and (^l]), to 
be —iJ2 P ^m, P { c n{t),c p (t)}. If we define the step func- 
tion u(t) as 



u{t) 



< 



1 , T > 

1/2 , r = 

, T < 



then the anticommutators for arbitrary r are given by 



{cn(t + r),F^(t)} = 2u(r) J2«p " iA^5 T , ){c n (* + r), c\{t)} 

p 

{c n (t + r), F^ t] \t)} = -M-r) J2^p + ^pSr.,o){c n (t + r), 4 (t)}. 

p 

I 



(55) 
(56) 



Finally, we show that Eq. ( |15| ) preserve the equal time 
anticommutators for the system operators. Since ( |l5| ) 
constitute a,n N x N system of equations, an explicit 
solution will be non-trivial for N > 2. However, it is 
already apparent that {c n (t), c m (t)} = since, as we 
stated before, c n (t) can be expressed in terms of the ini- 
tial states for the operators at to and {c n (to), c m (to)} = 
{c n (to), ak(to)} — 0. We are therefore left with showing 
that {c n (t) , c\ n {t)} = 5, hm . Calculating the derivative of 
the anticommutator using (J15|) and (p5|), one obtains, 



dt 



Integrating this expression and using the initial condi- 



{ Cn (t), <4(t)} = -i(n n - n m ){ Cn (t), cl(t)}. 



tion, {c n (t ), c^(to)} 
suit. 



we obtain the desired re- 



VIII. APPENDIX B: TUNNELLING COUPLING 
CONSTANTS 



In this appendix we derive explicit expressions for the 
tunnelling coupling constants K, n ^ and k n ^- Tunnelling 
in many-body systems was first treated by Bardeen 21] 
and later elaborated on by Prange |2(J and Harrison 2fj 
in the context of electron tunnelling across insulator junc- 
tions. Bardeen showed that if there exists a potential bar- 
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rier separating two regions of space, then the tunnelling 
matrix element, T a b, for the tunnelling of a particle from 
state 4>b(x) on one side of the barrier into the state <j> a {x) 
on the opposite side of the barrier is given by the off- 
diagonal current density in the barrier, 



2m 



h(x)—(t)* a (x) 



(57) 



Here x\ is a point inside the barrier and (fi a (x) and <fib{x) 
are the continuations of the single particle wave functions 
into the barrier where they decay exponentially. T a ,b is 
independent of the choice of x\ provided the energy dif- 
ference between the two states <fia and (fib is much less 
than the height of the barrier. The relationship between 
( p7|) and the overlap of the Hamiltonian between states 
localized on either side of the barrier is discussed in |20] . 

For the system illustrated in Fig. 1, we identify 
— ihn* n k with T a .b when <fi a {x) — (fi k (x) and (fib(x) = 
<fin\x). Similarly, —ihk^ k is equal to T a ^ when (f) a (x) — 

<j) k \x) and 4> b (x) = <fi { n\x). 

For the purpose of calculating the coupling constants, 
we take the <fin\x) to be the eigenstates of the finite 
potential well corresponding to d — > oo. In the classically 
allowed region, —a < x < a, <^4 (x) is proportional to 
cos(K n x) for n even and sin(K n x) for n odd while inside 

the barrier (x) ~ exp(— -\/2m(Vb — hn n )/h 2 \x\). The 
cavity wave numbers K n are determined by the solutions 
to 



K n a tan K n a = \/ (3 2 - (K n a) 2 . 
for n even and 



K n a cot K n a = - \J (5 2 - (K n a) 2 



(58) 



(59) 



for n odd @. 

For even n we find, 



with 

In a similar manner, we find for n odd, 



(61) 



with 



2K n 



iha 



m y L(l + (a/k) 2 ){K n a — tanif n a) 



e CTd sini^„a 



(62) 



Kn,k- (63) 

In both cases, a = ^/2mVo/h 2 ~ k 2 is the inverse of the 
penetration depth of the reservoir state into the barrier 
and 1 = 2ma 2 Vo / H 2 is the dimensionless barrier height. 

An interesting consequence of Eqs. (^lj) and (^) is 
that if n is even and m is odd or vice versa, then 



(64) 



iha 



2K n 



m V L(l + (a/k) 2 )(K n a + cot K n a) 



e ad cos K n a 



(60) 



On the other hand if n and m are both even or both odd, 
then 

(Tn,m "t" i^n^rn) ~t~ is)n,m ~1~ ^^n,m) — 2{^f n ra -\- i^ n m ). 

(65) 

Since (fin (x) is an eigenstate of the parity operator with 
parity (—1)™, it follows from Eq. ([?4j) that for a two- 
sided cavity only states of the same parity are coupled. 
This is a direct consequence of the our model system in 
Fig. 1 being invariant under spatial reflections, which 
implies that even when the coupling of the cavity states 
to the reservoirs is taken into account, parity will still be 
a good quantum number for the cavity states. No such 
result holds for the single-sided cavity since the cavity 
plus reservoir system is no longer invariant with respect 
to reflections. 

Finally we note that the approximation io k « 17„ used 
in deriving n ntk becomes exact for the evaluation of the 
7„ in and the cavity boundary conditions ( |27| ) and ( |3^ ) 
since in these expression n n ^ k is always multiplied by 
S(fl n - u k ). 
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